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Abstract 

A possible collective fermionic excitation in the ultrasoft energy-momentum 
region p < g 2 T is examined in Yukawa model with scalar coupling and quan- 
tum electrodynamics (QED) with g being coupling constant at extremely 
high temperature T where the fermion mass is negligible. We analytically 
sum up the ladder diagrams for the vertex correction in the leading order 
in QED, which is not necessary in the Yukawa model, and find that the 
fermion pole exists at u = ±|p|/3 — iC, with ultrasoft momentum p both for 
the Yukawa model and QED; £ is the sum of the damping rates of fermion 
and boson with hard momenta. We also obtain the expression of the residue 
of the pole, which is as small as of order g 2 . We show that the fermion prop- 
agator and the vertex function satisfy the Ward-Takahashi identity in QED. 
Thus we establish the existence of an ultrasoft fermionic mode at extremely 
high temperature, which was originally called phonino and was suggested 
in the context of supersymmetry and its breaking at finite T. We discuss 
the possible origin of such an ultrasoft fermionic mode without recourse to 
supersymmetry. The case of QCD is briefly mentioned. 



1. Introduction 

Revealing and clarifying the nature of possible quasiparticle and col- 
lective excitations is of basic importance for understanding a many-body 
system, in particular, in the low-energy regime. In fermion-boson systems 
such as Yukawa model, quantum electrodynamics (QED), and quantum 
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Figure 1: Schematic illustration of the energy hierarchy of bosonic and fermionic modes 
in a massless fermion-boson system at finite temperature: The vertical axis is the energy 
scale. The left and right part display the bosonic and fermionic modes, respectively. The 
bosonic modes may be of scalar-, gauge-bosonic and hydrodynamic ones. No established 
fermionic mode had been known in the ultrasoft region. 



chromodynamics (QCD) at so high temperature T that the masses of the 
particles are negligible, the average inter-particle distance is proportional 
to 1/T while some collective effects can be expected in the soft momentum 
scale 1/gT with g being the coupling constant. Indeed, soft bosonic modes 
in the longitudinal as well as the transverse channels exist and are known 
as plasmon [1], while the fermionic counter part is known as plasmino [2, 3], 
both of which have masses of order gT. In the lower energy region, there 
exist hydrodynamic modes of bosonic nature, which are actually the zero 
modes associated with the conservation of energy-momentum and charges. 
The energy hierarchy of bosonic and fermionic modes at high temperature 
may be summarized schematically as shown in Fig. 1. From the figure, one 
may have an intriguing but natural question whether such ultrasoft (< g 2 T) 
or zero modes can exist also in the fermion sector, possibly when the fermion 
system has a peculiar symmetry such as chiral symmetry. In this paper, we 
argue and demonstrate that such a fermionic mode can exist in this infrared 
energy region at extremely high temperature. 

Here we should mention that there have been some suggestive works for 
supporting the existence of such an ultrasoft fermionic mode at finite tem- 
perature. Historically, the ultrasoft fermionic mode at finite T was found 
in supersymmetric models as Nambu-Goldstone fermions called phonino as- 
sociated with spontaneously breaking of supersymmetry at T ^ 0, which 

2 



was shown by using Ward-Takahashi identity and a diagrammatic tech- 
nique [4, 5]. Here we note that the analysis in [4, 5] was performed in the 
temperature region much below the critical temperature, T c ~ m/g, where 
m is the bare mass. It implies that their analysis is only valid for gT -C m, 
but not for m <C gT for which our analysis in the present paper is con- 
cerned. The analysis was extended to QCD at so high temperature that the 
coupling constant is weak [6], in which a supersymmetry is still assigned at 
the vanishing coupling, and hence, the supersymmetry is, needless to say, 
explicitly broken by the interaction. Thus, there exists no exact fermionic 
zero mode but only a pseudo-phonino does. Although these analyses [4-6] 
are suggestive, it is still obscure whether a genuine ultrasoft fermionic mode 
exists when supersymmetry is absent, in particular, at extremely high tem- 
perature. 

Here we note that there have been suggestions of the existence of ul- 
trasoft fermionic mode at finite T even without supersymmetry. It was 
shown in one-loop calculations [7, 8] that when a fermion is coupled with 
a massive boson with mass m, the spectral function of the fermion gets to 
have a novel peak in the far-low-energy region in addition to the normal 
fermion and the plasmino, when T ~ m, irrespective of the type of boson; 
it means that the spectral function of the fermion has a three-peak structure 
in this temperature region. Recently, the present authors [9] have suggested 
that such a three-peak structure may persist even at the high temperature 
limit in the sense m/T — > 0, for the massive vector boson on the basis of 
a gauge-invariant formalism, again, at the one-loop order. Thus, one may 
expect that the novel excitation may exist in the far-infrared region also 
for a fermion coupled with a massless boson, although the one-loop anal- 
ysis admittedly may not be applicable at the ultrasoft momentum region. 
There are also the works suggesting the existence of the ultrasoft fermionic 
mode using the Schwinger- Dyson equation [10]; we note, however, that it is 
difficult to keep gauge symmetry in the Schwinger-Dyson approach at finite 
T. 

Now let us give a generic argument supporting the existence of an ul- 
trasoft fermionic mode at finite T on the basis of the symmetry of the self- 
energy for a massless fermion. In this case, the important point is that the 
fermion is chiral; the real part of the retarded fermion self-energy, T, r (uj, 0), 
at zero spatial momentum and the vanishing chemical potential, is an odd 
function of p° : 

Re£ fl (-p°,0) = -ReE K (p°,0). (1) 
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If the Re S R (p°, 0) is a smooth function at p° = 0, then, Re E R (— p°, 0) = 
at p° = 0, which implies ReG _1 (p°,0) = at p° = 0; the spectrum has 
a peak at the origin provided that the imaginary part of the fermion self- 
energy is not too large. This argument suggests that the existence of the 
ultrasoft pole may be a universal phenomenon at high temperature in the 
theory composed of massless fermion coupled with a boson. 

It is, however, not a simple task to establish that fermionic modes exist 
in the ultrasoft region on a general ground beyond the one-loop order ac- 
curacy because of the infrared divergence called pinch singularity [6, 11-14] 
that breaks a naive perturbation theory, as will be briefly reviewed in the 
next section. We remark that the same difficulty arises in the calculation of 
transport coefficients [12, 14] and the gluon self-energy [13] in the ultrasoft 
energy region. Therefore, in this paper, we analyze the fermion propagator 
in the ultrasoft energy region in Yukawa model and QED using a similar 
diagrammatic technique in Refs. [5, 6, 12, 14] to regularize the pinch sin- 
gularity. We shall show that the retarded fermion propagator has a pole 
at p° = ±|p|/3 — iC, (C ~ Q A T\iig~ l for Yukawa model and ~ g 2 T In g~ x 
for QED) with the residue Z ~ g 2 for ultrasoft momentum p taking into 
account the ladder summation. 

This paper is organized as follows: In Sec. 2, we discuss the ultrasoft 
fermionic mode in Yukawa model as a simple example without supersym- 
metry. In Sec. 3, we examine the ultrasoft fermionic mode in QED. We 
analytically sum up the ladder diagrams giving the vertex correction in the 
leading order, and find the existence of the ultrasoft fermionic mode as in 
the Yukawa model. We shall also show that the constructed propagator 
and the vertex satisfy the Ward-Takahashi identity. Section 4 is devoted to 
the discussion on the physical origin of the ultrasoft fermionic mode, and 
in Sec. 5, we will summarize our results. 

2. Ultrasoft Fermionic Modes in Yukawa model 

Let us start with Yukawa model, which is the simplest model to study 
the ultrasoft fermionic modes. Generalization to gauge theory will be dis- 
cussed in Sec. 3. As stated in Introduction, we deal with a massless fermion 
(denoted by ip) coupled with a scalar field through the interaction La- 
grangian Cj = gcf)^. We do not include the possible self-coupling of the 
scalar fields for simplicity. We calculate the fermion retarded self-energy 
and obtain the fermion retarded Green function with an ultrasoft momen- 
tum p < g 2 T. We first see that the naive perturbation theory breaks down 
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Figure 2: Diagrammatic representation of the fermion self-energy in Eq. (2) at one-loop 
order. The solid and dashed line are the propagator of the fermion and the scalar boson, 
respectively. In Eq. (14), the fermion and the boson propagator in the internal lines are 
replaced by the dressed propagators given in Eq. (10). 



in this case. Then, we shall show that a use of a dressed propagator gives 
a sensible result in the perturbation theory and that the resulting fermion 
propagator has a new pole in the ultrasoft region. 

2.1. One-loop calculation 

The retarded self-energy in the one-loop level is given by 

E&»(p) = *9 2 J-0jI [D S (-k)G*(p + k) + D*(-k)G s (p + k)] , (2) 

where D^ s {-k) and G^ s (p + k) are the bare propagators of the fermion 
and the scalar boson defined as 

Gg(k) = Q - n F (k°)\ ijk(27r)sgn(k°)6(k 2 ), (4) 

Do (*) =(\ + M* )) i(27r)sga(k°)5(k 2 ). (6) 

Here, n F (k°) = l/(exp(k°/T) + 1) and n B (k°) = l/(exp(k°/T) - 1) are the 
Fermi-Dirac and Bose-Einstein distribution functions, respectively. In the 
present analysis, we have preferentially employed the real-time formalism 
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in Keldysh basis [15]. The diagrammatic representation of Eq. (2) is shown 
in Fig. 2. Inserting Eqs. (3) through (6) into (2), we obtain 



J barc 



X p2 + 2 P -flto +p o )£ Wsgn(^ 
1 n F (k°+p^ * + * 



(7) 



2 v ' J p 2 + 2p-k + ik°e 

x (2vr)sgn(A; + p°)6((k + p) 2 )] , 

where we have used the on-shell conditions for the bare particles, k 2 = 0, 
and (k+p) 2 = in D(j(—k) and G(j(p+k). Then, for small p, the self-energy 
is reduced to 

74/ 



2p-k + ik°e ' 
where 

K(k) = (2tt) sgn^ )^ 2 )^^ ) + n B (A; )). (9) 

Note that K(k) is independent of p. This approximation is equivalent to 
the HTL approximation [2]. The HTL approximation is, however, only valid 
for p ~ gT, and not applicable in the ultrasoft momentum region. In fact, 
the retarded self-energy in the one-loop level obtained with use of the bare 
propagators is found to diverge when p — > 0, since the integrand contains 
1/p-k. This singularity is called "pinch singularity" [6, 11-14]. 

The origin of this singularity is traced back to the use of the bare prop- 
agators because the singularity is caused by the fact that the dispersion 
relations of the fermion and the boson are the same and the damping rates 
are zero in these propagators. For this reason, one may suspect that this 
singularity can be removed by adopting the dressed propagators taking into 
account the asymptotic masses and decay widths of the quasiparticles, as 
will be shown to be the case shortly. 

Since the leading contribution comes from the hard (k ~ T) internal and 
almost on-shell (k 2 m 0) momentum 1 , we are led to employ the following 



^^Here we note that the case where the internal momenta are soft (k ~ gT) or smaller 
is not relevant: In fact, the HTL-resummed propagators [16] should be used for soft 
momenta. However, the dispersion relations of the fermion and the boson obtained from 
these propagators arc different from each other, so the pinch singularity will not appear 
in this case. 
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dressed propagators for the fermion and boson: 



GR(k ^- v-4w (io) 

O\k)^{l-n F ^)) t {k2 _^ k ° iQ(k0r (11) 

(12) 



where rrif = gT /(2\/2) and m& = gT/y/6 are the asymptotic masses of 
the fermion and the boson at fc 2 ~ 0, respectively [17-19]. The damping 
rates of the hard particles, (f and (b, are of order g A Thxg~ x . The loga- 
rithmic enhancement for the damping rate is caused by the soft-fermion 
exchange, which is analogous to that of the hard photon [20]. Note that 
these resummed propagators are the same as those used in [6], except for 
the smallness of the damping rates: We remark that such a smallness is not 
the case in QED/QCD, where the damping rate is anomalously large and 
of order g 2 T\ng~ l ("anomalous damping") [21]. 
Using these dressed propagators, we obtain 

^ - 1 (0^> i + 2 ,W < 14 > 

for small p, where 5m 2 = m 2 . — m 2 = g 2 T 2 /24, £ = Q + K(k) = 
(g 2 1 8m 2 ) K (k) , and p M = (p° + i(,p). We have used the modified on-shell 
condition of the quasi-particles, k 2 —m 2 +2i(fk° = and k 2 —m 2 .+2i(fb° = 0, 
to obtain the denominator of the integrand in Eq. (14). We have also 
neglected rrib, rrif, (b, and (f in K(k), since the leading contribution comes 
from hard momenta k ~ T. It is worth emphasizing that thanks to 5m 2 
and (, T, R (p) given in Eq. (14) does not diverge in the infrared limit, p — > 0. 

Before evaluating Eq. (14), we introduce the the following dimensionless 
value: 

f d A k ~ n 2 T 2 

which is of order unity. This value will characterize the strength of residue 
of the pole for both Yukawa model and QED. 
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We expand the self-energy in terms of p M instead of p M itself. This is the 
key point of our expansion, which enable us to analytically find the pole of 
the ultrasoft fermionic mode. Then, the leading contribution is 

* - / w^ km2 i^ = 4 ((p0 + <)7 ° + vp ' 7) ' (16) 

with Z = (yf 2 /(8A 2 7r 2 ) and v = 1/3. Note that the zeroth-order term is 
absent, which implies that there is no mass term. Thus, we obtain the 
fermion propagator in the ultrasoft region as 

G R (p° p) = - - 

Z f 7 — p • 7 7° + p • 7 



2 \p° + v\p\ + iC p° — v\p\ + z£ 

Here we have decomposed the fermion propagator into the fermion number 
+1 and —1 sectors in the second line. These two sectors are symmetric 
under p -H- — p and v <H> —v, so we analyze only the fermion number +1 
sector in the following. 

From Eq. (17), we find a pole at 

p° = -v|p|-iC. (18) 

Note that the real part of p° is negative for the fermion sector, which sug- 
gests that this peak has an antifermion-hole-like character like the antiplas- 
mino [3]. The dispersion relation of the real part, Rep = — f|p|, is shown in 
the left panel of Fig. 3 together with the HTL results [2, 3] for comparison, 
where the coupling constant is chosen as g = 0.1. The imaginary part of 
the pole reads 

£ ~ g *T hi g- 1 , (19) 

which is much smaller than those of the normal fermion and the antiplas- 
mino [17]. Since the real part and the imaginary part of the pole are finite 
for |p| ^ 0, this mode is a damped oscillation mode. The residue of the 
pole is evaluated to be 

Z =8^=7^~^ < 20) 

which means that the mode has only a weak strength in comparison with 

those of the normal fermion and the antiplasmino, whose residues are order 
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\v\/{g 2 T) P °/(9 2 T) 

Figure 3: Left panel: The dispersion relation in the fermionic sector. In both of the two 
figures, the coupling constant is set to g = 0.1. The vertical axis is the energy p°, while 
the horizontal axis is the momentum |p|. The solid (blue) lines correspond to the normal 
fermion and the antiplasmino, while the bold solid (red) one to the ultrasoft mode. 
Note that since we focus on the fermion sector, the antiplasmino appears instead of the 
plasmino. The dotted lines denote the light cone. Since our analysis on the ultrasoft 
mode is valid only for |p| <C g 2 T, the plot for |p| > g 2 T may not have a physical meaning. 
The residue of the antiplasmino becomes exponentially small for |p| ^> gT, so the plot 
of the antiplasmino does not represent physical excitation for |p| ^> gT, either. Right 
panel: The spectral function in the fermion sector, Eq. (21), as a function of energy p° 
at zero momentum. 



of unity. It is worth mentioning that such smallness of the residue is actually 
compatible with the results in the HTL approximation: The sum of the 
residues of the normal fermion and the anti-plasmino modes obtained in the 
HTL approximation is unity and thus the sum rule of the spectral function 
of the fermion is satisfied in the leading order. Therefore, one could have 
anticipated that the residue of the ultrasoft mode can not be the order 
of unity but should be of higher order. Equations (17) through (20) for 
Yukawa model with a scalar coupling are obtained for the first time. 

The pole given by Eq. (18) gives rise to a new peak in the spectral 
function of the fermion as 

p + (p°, p) = -Im — ~\ (21) 
7T p u + v\p\ + ^C 

which is depicted in the right panel of Fig. 3, where |p| is set to zero. Since 
the expression of ( for the Yukawa model is not available in the literature, 
we simply adopt ( = g A T 'In g~ l / (2ir) in the figure. 
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2.2. Absence of vertex correction in Yukawa model 

So far, we have considered the one-loop diagram. We need to check 
that the higher-order loops are suppressed by the coupling constant. This 
task would not be straightforward because, 5m 2 ~ g 2 T 2 appears in the 
denominator, as seen in Eq. (14), which could make invalid the naive loop 
expansion. The possible diagrams contributing in the leading order are 
ladder diagrams shown in Fig. 4 because the fermion-boson pair of the 
propagators gives a contribution of order 1/g 2 , and the vertex gives g 2 . 
However, there is a special suppression mechanism in the present case with 
the scalar coupling. 

For example, let us evaluate the first diagram in Fig. 4, at small p. The 
self-energy is evaluated to be 

l K ( k ) / 77^ K ( l > ,nu n T~i • ( 22 ) 



(2vr) 4 v ' J (2tt) 4 w (2k- 1) 5m? 

Since there are four vertices and two pairs of the propagators whose mo- 
menta are almost the same, the formula would apparently yield the factor, 
K{k)K{l) ~ g A x (5m 2 )~ 2 ~ g°. One can easily verify that this order es- 
timate would remain the same in any higher-loop diagram, so any ladder 
diagram seems as if to contribute in the same leading order as explained. 
However, this is not the case for Yukawa model with the scalar coupling. 
An explicit evaluation of the numerators of the fermion propagators gives 
— f)f = t/k 2 — I/tl 2 , which turns out to be of order g 2 . This is because 
the internal line is almost on-shell, i.e., k 2 , I 2 ~ g 2 T 2 , which comes from the 
asymptotic masses squared. An analysis shows that the same suppression 
occurs in the higher-order diagrams such as the second diagram in Fig. 4. 
Thus, the ladder diagrams giving a vertex correction do not contribute in 
the leading order in the scalar coupling, and hence, the one-loop diagram 
in Fig. 2 with the dressed propagators solely suffices to give the self-energy 
in the leading-order. 

We remark that a similar suppression occurs in the effective three-point- 
vertex at p ~ gT [17]. We also note that this suppression mechanism is 
quite similar to that found in a supersymmetric model for an intermediate 
temperature region in the sense that gT <C m [5], whereas we are dealing 
with extremely high-T case. It should be emphasized that this suppression 
of the vertex correction is not the case in QED/QCD, where all the ladder 
diagrams contribute in the leading order and must be summed over [6], as 
will be shown for QED in the next section. 
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Figure 4: Some of the ladder class diagrams. The solid and dashed line are the dressed 
propagators of the fermion and the scalar boson, respectively. 



3. Ultrasoft Fermionic Modes in Quantum Electrodynamics 

Next we explore whether the ultrasoft fermionic mode also exists in 
QED at high T. One might expect that the analysis would be done in much 
the same way as in the Yukawa model. It turns out, however, that the 
analysis is more complicated and involved. It is necessary to sum up the 
contributions from all the ladder diagrams even apart from the complicated 
helicity structure of the photon. In this section, we successfully perform 
the summation of the ladder diagrams in an analytic way, and obtain the 
fermion propagator that is valid in the ultrasoft region. Then we evaluate 
the pole in the ultrasoft region explicitly and examine the properties of the 
ultrasoft fermionic mode in QED. We also discuss whether the resummed 
vertex satisfies the Ward-Takahashi identity. 

3.1. One-loop calculation 

First, we evaluate the contribution from the one-loop diagram. We use 
g as a coupling constant of QED instead of the standard notation e to make 
it clear that the same order counting appears as that in the Yukawa model. 

The dressed propagators with hard momenta read 



G R (k) 



-ft 



(23) 



Dt{k) 



k 2 - m 2 e + 2iC e fc°' 



(24) 



P-m2 + 2< 7 fc ' 




(25) 



(26) 
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Here m e = gT/2 and m 7 = gT/\^6 are the asymptotic masses of electron 
and photon respectively [18, 19]. The damping rates of electron and photon 
are estimated as ( e ~ g 2 T In g~ x [21] and C 7 ~ g^Tlng^ 1 [20]. Note that ( e 
is much larger than that in the Yukawa model, which is called "anomalous 
damping" [21]. This large electron damping makes the damping rate of the 
ultrasoft mode larger than in the Yukawa model. P^ u {k) is the projection 
operator on the transverse direction, 

p ^( k ) = 9tii9uj(Sij - kkj), (27) 

with k = fc*/|k|. Here we have adopted the Coulomb gauge, in which the 
analysis becomes simple thanks to the transversality of the photon propaga- 
tor. We note that the u^u v /|k| 2 term, where w M = (1, 0), has been omitted 
in Eqs. (24) and (26) because that term vanishes after the k° integral. The 
analysis in other gauge-fixing conditions are more complicated, and will be 
reported elsewhere [22]. 

By using these resummed propagators, the one-loop contribution in the 
ultrasoft region is evaluated as 

^LioopW = i9 2 j -0^1" [Dl{-k)G R {p + k) + D^(-k)G s (p + k)]Y 

d * k K(k) £ 
(2tt) 4 1 ] l + 2p- k/Sm 2 ' 

(28) 

where Sm 2 = m 2 — ml, ( = Ce + C7 — Ce- Here we have used the same 
notation for 5m 2 and ( as those in the Yukawa model, although their para- 
metrical expressions are different from each other. The factor two in the last 
line of Eq. (28) comes from two degrees of freedom of photon polarization. 
At the one-loop order, we obtain 

^Lioc>) = ~ 1 ^ 2 ((P° + iC)7° + • 7) • (29) 

We note that this expression has the same structure as that for the Yukawa 
model; see Eq. (16). 

3.2. Ladder summation 

As already mentioned, the ladder diagrams contribute to the leading- 
order unlike in the Yukawa model. In this subsection, we sum up all the 
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ladder diagrams, and obtain the analytical expressions of the pole position 
and the residue of the ultrasoft mode. 

For this purpose, we introduce the vertex gT^^p, k) defined through the 
following self-energy: 

h W " J (2tt)< W l + ip-k/Sw? ■ (M> 

Here the vertex contains the contributions from all the ladder diagrams and 
satisfies the following self-consistent equation: 

dAk ' r>nn v + $ + + 



(p + k + k') 2 l + 2p- k' /5m 



2 



xpT(k')re(p,k'). (31) 

Equations (30) and (31) are represented diagrammatically in Fig. 5. We 
have used the same approximation as that used in the derivation of Eq. (8) 
for the propagator of fermion and photon. We should remark here that this 
summation scheme using the self-consistent equation was first constructed 
in Ref. [6]. However, we also note that the equation has never been solved 
either analytically nor numerically. In the following, we solve this self- 
consistent equation analytically for small p, and show that the dispersion 
relation does not change from that in the one-loop order even after incor- 
porating all the ladder diagrams. 
At small p, Eq. (31) reduces to 

rn k) _ f d*k> Pj p {k>) + k^nlr^p, k>) 

l IP, A;) -7 +J i 2n)* K{k) l + 2(p.k<)/5m* k ■ k> 

(32) 
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where we have dropped the term proportional to ^ because it only gives 
higher order contribution after being multiplied by f, in the numerator of 
Eq. (30), because f = k 2 ~ g 2 T 2 . 

Let us solve the self-consistent equation (32). We expand the vertex 
function as 

r»(p,k) = r£(k) + 6r»(p,k), (33) 

where r^{k) is of order unity and 5r^{p, k) is of order p/(g 2 T). 
We first evaluate r^(k), which can be decomposed as follows: 

rg(k) = YMk) + k^B(k) + u"C{k), (34) 

where A, B, and C are 4x4 matrices. Then the self-consistent equation 
for r^k) becomes 

/d A k' ~ k'k v P T (k' s \ f y p 
j^yA k ') kk' ^ 

(35) 

where B(k) and C(k) in the right hand side vanish due to transversality of 
the photon propagator: P^ u {k)k u = P^ u {k)u v = 0. By assuming that A(k) 
is a constant, the integral becomes 



(2tt) 4 v ' 1 k-k' V k ° 

where we have dropped f, term and imposed k 2 = as before. Then, we 
find 

A=—l, m = ¥TTt f, C(*) = 0, (37) 

where 1 is the unit 4x4 matrix. 

Next we evaluate 5T(p, k). From Eq. (32), expanding -T M (p, k) in terms 
of 2p ■ k' /5m 2 , we find the self-consistent equation for 5T(p, k) as 

(38) 



x 



- A y^Jr + 5rP iv, k>) 
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It is not easy to analytically solve Eq. (38) directly. So we instead calculate 
the following function: 

Sn(p) = j ^K(k)P^(k)r^(p,k). (39) 
Then Eq. (38) leads to the following closed equation, 



k ■ k' 5m 2 



+ 



d A k' 



ri 4 k ~ k' v 

" K{k)P T (k)rh 



(2tt) 4 v ' ^ ' 1 ' k ■ k 

xk(k')pJ x (k')Yf6r\ P ,k f ) 
(p)-xsn(p). 



(40) 



Here we have used Eqs. (36) and (39) in the second line. The solution to this 
equation is readily found to be SII(p) = — \A 2 U onc _\ oop (p) , where Eq. (37) 
is used for A. Then, the self-energy is evaluated as 

= ^one-loop(p) + Sn{p) (41) 

= -^(7V + »C) + «P > 7), 



where the residue is 

f 

16vr 2 A 2 



Z = t/tt,(1 + A) 2 . (42) 



The pole of the ultrasoft fermionic mode has the velocity v = 1/3, damping 
rate (, and the residue Z. The dispersion of the mode is the same as that in 
the one-loop level whereas the residue is changed. This is our main result 
for QED. 

3.3. Ward-Takahashi identity 

In this subsection, we examine whether the summation scheme, Eqs. (30) 
and (31), is consistent with the U{1) gauge symmetry. Concretely, we 
check that our resummed vertex function and self-energy satisfy the Ward- 
Takahashi (WT) identity in the leading order of the coupling constant. 
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The WT identity reads 



k^r^p, k) = i> + $ - E R (p + k)-p* + E R (p). (43) 

Since r^p, k) contains two separated scales k ~ T and p < g 2 T, we need 
to treat them carefully. For the hard part, E R (p + k) ~ E R (k) is of order 
g 2 T, which is negligible compared to ^. In addition, the momentum depen- 
dent part E R {jp + k) — E R {k) ~ g 2 p is also negligible compared to E R {p). 
Therefore the WT identity reduces to 

k^(p,k) = k' + Z R (p). (44) 

On the other hand, multiplying Eq. (31) by k^, we have 

d A k' ~ v t + ft + t' 



k ll r»( P ,k) = u + t)-t+ J 



(p + k + k') 



2 



{{i>+n^)-i>-m+f) pT{k r )rP{v k>] 

X 1 + 2p • Jfe'/Jm* ^ {k}1 {P ' k} 



^ + r R ( P ), 



(45) 



where we have dropped the terms of order g 2 T, and used Eq. (30) in the last 
line. This expression coincides with Eq. (44), so our self-consistent equation 
satisfies the WT identity. We note that this proof was made without using 
the expansion in terms of p/g 2 T. For this reason, Eq. (44) is generally valid 
for p < g 2 T, as well as for p <C g 2 T. 

Next, we check whether the explicit solution Eq. (37) of the self-consistent 
equation (31) at zeroth order in p/(g 2 T) satisfies the WT identity to see 
the consistency with the gauge symmetry of the following two conditions 
adopted to obtain Eq. (37): One is that terms proportional to ft in r(p, k) is 
dropped, since they are negligible in the self-energy due to ft ft = k 2 ~ g 2 T 2 . 
The other is that we imposed the on-shell condition k 2 = 0, because the 
internal photon in the self-energy is almost on-shell. Using the same con- 
ditions, we expect that the vertex function satisfies the WT identity in the 
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leading order of the coupling constant. In fact, we have for the zeroth order 
in p 



k^{k) ~ ftA + k 2 B(k) ~ 0. (46) 

Here we have dropped f, and k 2 in the last equality. 

Finally, we check that the equation determining the vertex function 
Eq. (38), which is first order in p, satisfies the WT identity. By multiplying 
this equation by fc M , we obtain 



(47) 



where we have dropped ^ as in the previous equation and used Eq. (41). 
Therefore, our analytic solution of the self-consistent equation satisfies the 
WT identity in the leading order of the coupling constant. 

We note that without the summation of the ladder diagrams, the WT 
identity is not satisfied when the external momentum of fermion is ultrasoft. 
By contrast, the ladder summation was unnecessary in the Yukawa model, 
in which the gauge symmetry is absent. 



4. Discussion on Origin of Ultrasoft Fermionic Mode 

Let us discuss the physical origin of the ultrasoft mode for clarifying its 
possible universality or robustness at high temperature. 

As mentioned in Introduction, an ultrasoft fermionic mode with a van- 
ishing mass can appear as a phonino associated with the spontaneously bro- 
ken supersymmetry at finite T [4-6], below T c : In the massive Wess-Zumino 
model, the supersymmetric cancellation of the fermion mass is essential to 
make the phonino excitation [5]. In fact, the massless fermionic mode is 
realized as a consequence of an exact cancellation of the self-energy with 
the finite bare mass at the vanishing external momentum due to the spon- 
taneous breaking of supersymmetry. This is quite in contrast to the case of 
the Yukawa model dealt in the present work, where the ultrasoft fermionic 
modes do not appear when the fermion mass is large. This means that 
the ultrasoft fermionic mode like the phonino can appear even in a non- 
supersymmetric model like the Yukawa model, and there the masslessness 
of the fermion is an essential ingredient in realizing it. Nevertheless it is 
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interesting that such an exotic fermionic mode can appear both in super- 
symmetric and nonsupersymmetric models with quite different mechanisms 
at high temperature. 

For getting further intuition into the possible mechanism for realizing 
the ultrasoft fermionic mode in the Yukawa model 2 , we note that Kitazawa 
et al. argued that the level repulsion due to the Landau damping causes the 
three-peak structure in Ref . [7] , though in the case where a massless fermion 
is coupled with a massive boson and thus the fermion-boson mass difference 
squared 5m 2 is nonzero. Indeed a finite Sm 2 plays an important role in the 
appearance of the ultrasoft mode even when the boson is massless because 
it ensures the smoothness of the self-energy at the origin: As discussed in 
Introduction, the real part of the self-energy vanishes at the origin from the 
symmetry for the massless fermion. The residue of the pole is proportional 
to (Sm 2 ) 2 and vanishes if m/ = m& in the Yukawa model or m e = m 7 in 
QED. In the present work, the difference in the masses squared becomes 
finite because the effects beyond the HTL approximation is taken into ac- 
count and found to be 5m 2 ~ g 2 T 2 . The mass difference in turn causes the 
smallness of the imaginary part: The imaginary part of the self-energy can 
originate from the boson emission, the Landau damping, and the imaginary 
part of the dressed propagators obtained beyond the HTL approximation. 
The former two contributions are found to be zero at the origin [7, 9] due to 
the nonzero mass difference or kinematics at the one-loop order. Therefore, 
the leading contribution of the imaginary part is solely given by the damp- 
ing rates of the hard particles. Thus the ultrasoft fermionic mode gets to 
exist without supersymmetry at high temperature. 

5. Summary and Concluding Remarks 

We have investigated the spectral properties of massless fermion with 
the ultrasoft momentum (p < g 2 T) in Yukawa model and QED in the 
Coulomb gauge at high temperature. We have first indicated that a re- 
summed perturbation theory [5, 6] is needed beyond the conventional HTL 
approximation to get a sensible spectral properties in the ultrasoft region: 
For a consistent calculation in the relevant order of the coupling, we have 
shown that the use of the dressed propagators with the asymptotic masses 



2 The mechanism of realizing the ultrasoft fermionic mode in QED is expected to the 
same as that in the Yukawa model because the multi-loop diagrams in QED affect only 
the expression of the residue. 
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Tabic 1: The expressions of the dispersion relation, the damping rate, and the residue of 
the ultrasoft mode in the Yukawa model and QED. 



Yukawa model 



QED 



dispersion relation 
damping rate 
residue 



±IpI/3 
Cf + Cb ~ 9 A T In g 



-i 




±IpI/3 
- g 2 Tlng 
2 /(144tt 2 ) 



and the decay widths for the both fermion and boson is necessary, and the 
vertex corrections due to the summation of ladder diagrams is to contribute 
in the leading order. In QED, they turn out to be of leading order and are 
non-negligible, while it is not the case in the Yukawa model. Consequently, 
the WT identity is satisfied in QED. 

Thereby we have established that the resulting fermion propagator de- 
velops a pole at ultrasoft momentum, which was first suggested in [5] for 
supersymmetric models in different temperature region from ours, and in [6] 
for gauge theory in the same temperature region as ours: Its pole position, 
width, and the residue have been obtained for the first time in the present 
work as a function of the asymptotic masses and damping rates; they are 
summarized in Table. 1. We note that we have obtained the pole of the 
propagator by expanding the propagator around p° = |p| = instead 
of p° = |p| = in QED, where the damping rate of the hard electron is 
anomalously large [21] unlike the Yukawa model. If one expands the prop- 
agator around p° = |p| = 0, it is assumed that p <C 5m 2 /T. The pole that 
we have found has the imaginary part Imp = ( ~ g 2 T, so the expansion 
would break down because pT/5m 2 ~ 1 at the pole position. We also re- 
mark that although the present analysis is given for a massless fermion, the 
ultrasoft fermionic excitation should exist even for a massive fermion, if the 
bare mass is smaller than g 2 T, which is the smallest scale in our analysis. 

We should also emphasize that the existence of a peak in the fermion 
spectral function at the origin implies that the fermion spectral function has 
a three-peak structure, although the central peak has only a small strength. 
The development of such a three-peak structure at intermediate tempera- 
ture was suggested in the case where the boson is massive irrespective of 
the type of the boson [7-9], as mentioned in Introduction. Thus, we see 
that the three-peak structure of the fermion spectral function is a robust 
phenomenon to be seen at both of intermediate and high T, at least in the 
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Yukawa model. In QED, further analysis is needed to establish the existence 
of the peak in the spectral function at high T. In general, the propagator 
can be expanded by the partial fractions [23], and several poles near the 
real axis contributes to the spectral function, although infinite number of 
poles generally exist in complex energy plane. The pole that we have found 
in this paper is one of them. It is, therefore, needed to examine whether 
the contributions from other poles are small 3 . 

Although the present analysis is restricted to the Yukawa model with a 
scalar coupling and QED, a similar analysis can be performed in QCD [6, 
22]. We note, however, that the self-coupling between the gluons leads to 
additional ladder diagrams that should be also summed up [6]. The detailed 
analysis will be reported elsewhere [22]. 

It is known that the HTL approximation is equivalent to the Vlasov 
equation [24]. For the bosonic channel in ultrasoft or softer momentum 
region, the ladder summation is required and corresponding kinetic theory 
is the Boltzmann equation [13]. It is interesting to develop a kinetic equa- 
tion for fermionic channel in the ultrasoft momentum region, which will be 
discussed in our future work [22]. 
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